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<D Abstract 

In the absence of errors, the dynamics of a spin chain, with a suitably engi- 
neered local Hamiltonian, allow the perfect, coherent transfer of a quantum 
state over large distances. Here, we propose encoding and decoding proce- 
dures to recover perfectly from low rates of systematic errors. The encoding 
and decoding regions, located at opposite ends of the chain, are small com- 
pared to the length of the chain, growing linearly with the size of the error. 
We also describe how these errors can be identified, again by only acting on 
qh the encoding and decoding regions. 

1 Introduction 

00 

£ — ■ Quantum state transfer, i.e., the coherent transfer of an unknown quantum 

state between two distant locations, plays a fundamental role in quantum 
information processing, facilitating the production of entanglement between 
^2 separated parties, and generating the long-range two-qubit interactions re- 

{NJ quired in a quantum computer. If the available inter-qubit interactions only 

act locally, one of the most effective ways to achieve this state transfer is to 
t> use a spin chain as a quantum wire [1, 2, 3] (see [4] for a review). In this 

scheme, the dynamics of a suitably designed Hamiltonian perfectly trans- 
!— j fer a quantum state between the two ends of an arbitrarily long chain. Its 

main strength is that external control is not required during the transfer 
while confining the encoding and decoding areas to the boundary regions of 
the chain, so that it is far superior to the naive solution of applying a se- 
quence of SWAP gates. Since the protocol is designed to operate in the same 
solid-state technology that the rest of the quantum device is fabricated in, 
interfacing with other technologies, such as using photons as 'flying qubits' 
[5], is unnecessary. 
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One of the promising features of a spin chain is that, in only requiring in- 
teractions with the ends of a chain, the majority of the spins can be shielded 
from the effects of noise. Nevertheless, errors will inevitably be generated 
in a practical implementation by manufacturing imprecisions and through 
coupling with the environment. The main effect of these errors is to compro- 
mise the transfer by introducing unwanted excitations and destroying the 
coherence. In order to reduce this effect, one can consider encoding in an 
error correcting code. However, applying standard error-correcting codes 
to state transfer is very inefficient; any error occurring during the transfer, 
even on a single qubit, is spread over the whole chain by the action of the 
Hamiltonian. Hence, one expects that a generic error-correcting code would 
require a large number of encoding qubits per logical qubit, conflicting with 
the central tenet of the state transfer protocol. It is therefore necessary to 
design ad hoc error correcting codes. In addition to the importance of error 
correction on the spin chain itself, such studies comprise the essential ini- 
tial studies of equivalent questions in systems where Hamiltonian dynamics 
are used in an information processing capacity, such as in a computational 
architecture [6]. 

To date, studies of error correction in state transfer have adopted a 
cause-based approach. Manufacturing errors have been modelled as random 
static or dynamical defects, [7, 8], or studied via simulations [9], and can 
be circumvented (at the cost of a significantly increased transfer time) by 
encoding across multiple parallel chains [10, 11]. The effects of timing er- 
rors can be minimised by optimally tuning the couplings between the chain 
qubits [ ]. The studies of environmental induced noise models are, to date, 
extremely limited, being based around the very crude attempts such as in 
[4, 13], or with very specific models with little physical motivation, [14]. 
Nevertheless, there are specific occasions when perfect state transfer can be 
recovered, albeit with a doubled transfer time [15]. No general procedure to 
protect against noise is yet known; indeed, aside from some of the arguments 
of [4], all these prior results have made use of a deliberately induced sub- 
space structure in the Hamiltonian which is assumed to be preserved, even 
in the presence of these errors 1 . One promising avenue, recently proposed 
in [16], advocates the use of topologically protected subspaces in chiral spin 
liquids in order to protect against decoherence during the state transfer. 

This cause-based classification has so far concealed an interesting cat- 
egory of errors, systematic errors, by which we mean errors that in each 

1 This is usually, but not exclusively, the spin preserving property that we describe in 
Sec. 1.1. 
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transfer experiment are drawn from a fixed set of possible errors (whose 
description includes position, time, and type of error). This key feature en- 
ables the detection, once and for all, of the set of possible errors, potentially 
allowing us to protect against them. In this paper, we develop a strategy to 
recover perfectly from low rates of systematic error on a spin chain, encod- 
ing the information in a protected set of states and decoding it accordingly. 
The number of qubits used for the encoding and decoding scales linearly 
with the sizes of the errors in the set. We start, in Sec. 1.1, by introducing 
the requisite information about spin chains, particularly concentrating on 
the mapping to a system of free fermions, the mathematical structure which 
we rely on throughout the paper. The encoding and decoding strategies 
are presented in Section 3 under the assumption that the errors are known, 
before we describe an example of how the errors can be identified in Section 
4, while also indicating that significant efficiency savings can be made over 
and above the sufficient conditions previously derived. 

1.1 Spin Chains and Fermions 

In an ideal state transfer protocol, we are provided with an unknown quan- 
tum state 

|X) = (a|0)+/?|l», 

and tasked with transmitting this state perfectly to a distant recipient, using 
only a fixed Hamiltonian with local couplings. The most commonly studied 
system is that of a spin chain, initially because this is the optimal geometry 
for transferring a state as far as possible. However, the many advantages of 
this choice include the ability to prove necessary and sufficient conditions for 
achieving perfect state transfer [4], and the ability to perform the Jordan- 
Wigner transformation on the system [J 7]. We will also take advantage of 
these benefits. 

In the simplest scenario, this state transfer is achieved by placing the 
state onto the first qubit of an N qubit chain, and initialising the rest of the 
system in a standard state such as lO}®^ -1 . Thus, the overall starting state 
is 

m = \x) ® io)^- 1 

By suitably selecting the properties of the Hamiltonian H, we evolve for 
some time, tf, known as the transfer time, and the corresponding unitary 
U(tf) = e~ lHt f creates the final output 

\^ F )=U(t f )\^ I ) = \0f N - 1 \x), 
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i.e. the state arrives on qubit N of the chain. We consider Hamiltonians of 
the class 

N-l N 
H = - Ji{XiX i+1 + YiY i+ x) -Y^BiZi 

i=l i=l 

where Xi,Y$ and Zj are the Pauli operators acting on the i-th qubit. This 
Hamiltonian has a very convenient subspace structure, described by 



N 



i=l 



0, 



meaning that the total number of ones (excitations) relative to the state 

|0) = |O) 0Ar 

is preserved. In particular, the one-excitation subspace, spanned by the 
states 

\n) = lO)®^- 1 ) |l)|0)® (JV - n) , 

is preserved. This allows us to describe the dynamics of this subspace by 
a Hamiltonian Hi, where (i?i)n,m = (r\ H\m). Perfect transfer is achieved 
by an appropriate choice of coefficients Jj and Bi, the necessary and suffi- 
cient conditions for which are described in [4], which result in the following 
properties: 

U(t f )\0) = |0) 

U(t f )\n) = exp(i(/>) \ N -n + 1 ) Vn = 1 . . . N 

for some real time tf and phase cf>. There are many such choices of coupling 
that achieve perfect state transfer. While the original [2, 3] has been shown 
to be optimal against a variety of criteria [12], we remain agnostic to the 
specific choice, and do not require more than the basic properties outlined 
here. The phase eft shall be neglected in the following, as it can be corrected 
by a local rotation. 

As we will be concerned by the action of errors which, generically, may 
introduce additional excitations during the state transfer, we shall switch 
to the fermionic picture, which lends itself to the description of variable 
excitation number. In this picture, the Hamiltonian reads: 

N-l N 
H = y~] J n (al t a n+ i + a) n+l a n ) + 2 ^ B n a} n a n 

n=l n=l 
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where 

This hopping Hamiltonian describes N non-interacting spinless fermions [4] . 

j n— 1 

The Jordan- Wigner transformation a\ = - Z m {X n — iY n ) , Z n = 1 — 

m=l 

2a\ia n provides the relation to the spin picture. Setting 0(t) = U{t)OU{t)\ 
for any operator O, the bilinear structure of H implies that 

N 

where (3 n ,m{t) = ex p( — iHit) \ m) [18] (i.e., the dynamics are entirely 
determined by the single excitation subspace, described by Hi). Hence, two 
fermions are perfectly transferred independently of one another, modulo an 
ordering phase: 

U(t f )a n a m UHt f ) 



1.2 Initial State 

We specified that for the basic state transfer protocol, the initial state of 
every qubit except the first one should be |0). However, a relatively mi- 
nor adaptation of the protocol makes this requirement unnecessary [13, 19]. 
Consider, instead, fixing the first two qubits to convey the information of 
the unknown state via an encoding 

(aa\ + (5a\) |00) , 

and the state of the rest of the chain can be completely arbitrary. The state 
of the rest of the chain is described by an arbitrary state 



where a x = 11^=3 ( a \) Xi uses the bits of the binary string x to specify which 
of the N — 2 qubits have creation operators acting on them. Hence, 



m = (aa\ + /3a\) \00) i V 7*4 |0) 




m-2 
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but by the properties of perfect state transfer and the non-interacting fermion 
picture, we immediately have that the output state is 

m = I £ 7**(-ir* + ("")4 ) |0)® JV - 2 («a] v + ^a] v _ 1 )|00), 

where x R is the bit string x reversed and w x is the Hamming weight of x, 
i.e. the number of creation operators that act in at- Clearly, this state is 
separable about the partition (1 . . . N— 2) and (N — 1, N), so we can recover 
the transferred state from the other end of the chain. By linearity, any 
mixture of these initial states can also be tolerated. 



2 Modelling systematic errors 

We represent a systematic error by a set of possible errors {£{\ acting at 
times {te}. Since the {aj} form a complete basis, each of these Si can be 
expanded as 

s rrii rii 
i=l j=l 3 j=mi+l 

where k^p are the different positions that the creation/annihilation op- 
erators act on. The {^{\ are arbitrary coefficients. The overall error £ is 
taken to be trace preserving. For simplicity of exposition, we shall work with 
just one error £ acting at time t, the effect of multiple errors being easily 
incorporated. The error is assumed to act instantaneously, or, at least, for a 
time much smaller than tf, so that otherwise the time-evolution is governed 
by U 2 . Hence the state at time tf is \iJ>f) = £(tf — t)U(tf) \tpi) , with 

-*)= j2 ^n a i« (*/-*) n %f (*/-*)■ ( 2 ) 

i=l j=l 3 j=m.i+l 

The fermionic operators a,i(tf — i)'s represent the (possibly non-orthogonal) 
fermionic modes affected by the error. 

Even without knowing the specific form of the error, the encoding pre- 
sented in Sec. 1.2 already tolerates a vast rage of errors; those for which all 
the errors £ acting at time t satisfying £ (tf — t) = Ylx R a ^ r ( or an y mixture 
thereof), if |0) is used as initial state, since this is entirely equivalent to the 

2 We comment on the relaxation of this assumption in the conclusions. 
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^ Encoding 




Figure 1: General schematic for the operation of a perfect state transfer 
system in the presence of a systematic error at time t via an initial encoding 
and corresponding decoding in small regions at either end of the chain. 

case of the different initial state. Our task is now to describe how to deal 
with errors that do not satisfy this property, under the assumption that we 
know what these errors are. 

3 How to counteract systematic errors 

In order to counteract these systematic errors, we shall encode the informa- 
tion about the state \x) in the first D qubits at one end of the chain, which 
we call the encoding region. The encoding will be defined so that perfect 
state transfer can be achieved by applying a unitary decoding operator Up 
at time tf, acting just over a decoding region of size D at the opposite end 
of the chain (see Fig. 1). We refer to the region outside the decoding region 
as the complement region, which has size D = N — D. The size D will be 
determined by the number of errors that we have to encode against, but 
should be considered small compared to the transfer distance, N. 
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3.1 State Encoding 

We propose to encode the information in the state 

m = (aQl + PQ{)\fh) (3) 
where the encoding operators 

D 

Ql = E( £ >! + tf°*)> a G {0, 1} , e? , r,? G C (4) 
i=i 

are to be determined to ensure that the decoding operation Ud exists; the 
state |Y>o) fulfils: 

Q a |Vo) = 0, Vae {0, 1} . (5) 

We require the {Q a }'s to represent two orthogonal fermionic modes, which 
imposes the conditions 

i=l 

t D 

=0 ^ 2(€?7 / ?'+7/?€?') = 0.Va€{0 ) 1} (7) 

i=l 

This ensures that the logical 0, Q\ \ipo), an d the logical 1, Qj |"0o)j ar e or- 
thogonal, normalised states, (V>ol QaQ^i \ipo) = $„„'■ The idea behind the 
proposed definition is that if the encoding operators correspond to fermionic 
modes different from the ones affected by the error, the orthogonality of 
the encoding states would be unaffected by the error, and so would the 
encoded information. In addition, since the Hamiltonian describes indepen- 
dent fermions, the transfer of the information-carrying excitations would be 
independent of the excitations introduced by the error, making it possible 
to recover the information by applying an appropriate decoding operation 
at time tf. We shall now formalise this intuition. 

3.2 State Decoding 

Let us rewrite \iPf), highlighting the action of the time-evolved error and 
encoding operators on the decoding region. To do this, we describe the 
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action of any given fermionic operator as a product of the part that acts on 
the decoding region, and the part that acts on its complement. 

at = f] ® 1® D , } 3 = X -Z^- x ®{X- iY) I® 5 -- 7 , Vj < D 

a] = Z® D ® /t , /t = l -zi- b - x ® (X - iY) ® Vj > D , 

with = Sij, {fi,fj} = (and similarly for the f/s). Hence, a 

fermionic operator evolved in time from when it acts up to the decoding 
time, tf, is expressed as 

a n (t) = F n (t) 1® D + Z® D ® F n (i) , 

D N 
with F n (t) = P n ,m{t -tf)f m , F n (t) = /3 n ,m(* -tf)f m .As we saw 

m=l m=_D+l 

in Sec. 2, the -F n (t) are essentially irrelevant, and it is the F n (t), the parts 
acting on the decoding region, that we must take care of. Substituting in 
Eq. (2), one obtains 

s 

£(tf-t) = J2li £ ^°A.®PW (8) 

i=l x6{0,l}"i 

summing over the string rr in order to convey within the Fx all the different 
combinations of F k i acting on the decoding region 

p * ] = n (4) n 1 > a * = n^)*'' 

j=l ^ ' j'=mi+l j=l 
with (obtained by substituting F fe i for Ft.) conveying the equivalent 

3 3 

information for the complement region 3 . The complement of x is denoted 
by x. e x conveys the parity of the ordering of the string x, as arises from 
imposing a standard ordering of the fermionic operators. Rewriting the 
encoding operators in the same formalism, one has: 

U(t F ) = (aQ (t f ) + Qi(t f )) |Vw> , |V „ t ) = U{t f ) |Vo> (9) 

where 

D 

Qi(t f ) = Z® D ®q a , ql = ^ + ^"/iV-i+l • 

i=l 

3 The argument t has been omitted for simplicity. 
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To complete the rewriting, we note that it is always possible to write \ipout) = 
| <f>i) <S> \~4>i), where \<pi) is any (not normalised) state defined over the 
complement of the decoding region, and, via Eq. (5), is a solution to 



,o = 0, 1 



(10) 



One can find this state as an eigenstate of the operator q\q$q\q\ with eigen- 
value 1. Provided D > 2, such eigenstates always exist (the size of the 
Hilbert space spanned by the solutions of this equation is 2 D ~ 2 ; hence |Vw) 
can be chosen in a Hilbert space of dimension 2 N ~ 2 , i.e., condition (5) is not 
too restrictive). In what follows, we shall work with \ip ou t) = |^} ® \ip) for 
the sake of simplicity, and more general states follow by linearity. 
The state at time tf can be rewritten as: 



Mf) 




4>l) <g> (a 



+ P 



i(0 



(11) 



with 



(i) 



P^qo I 



,(<) 

*-x 



"a: <7l 



6« A 



and U|) 

We shall now prove that if Eqns. (6), (7) and (10) hold, it is sufficient to 
impose that 



t F t 

q a i r k (i) 



D 



(12) 



for both encoded states (a = 0, 1), and each of the possible operators Fj (j = 
1 . . . m) from all possible error strings (i = 1 . . . s), to ensure the existence of 
least one unitary Up that, applied at time tf, perfectly recovers the encoded 
information. Fulfilling these conditions then defines the coefficients of the 
encoding operators in Eq. (4), completely specifying our strategy. 



-. x e {0,l} n > Vi = 1. 



Let us introduce the sets of vectors: Z a = j 

a = 0, 1, representing the domain of £7d . Conditions (12) imply the crucial 
property 



P (i) 



la 



■\Y*q a P® , P^qt = (-l) w *qM ] , Va , V* , Vi . (13) 



Consequently, via (7) and (10), no matter what the error, the two logically 
encoded states remain unambiguously distinguishable, (Ox \ly ) = , Vx G 
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{0, l} ni , \/y 6 {0, l} n i , Vi, j. This is already sufficient to imply the existence 
of Up, although we give an explicit construction in the Appendix. 

We have just proven that conditions (6), (7), (10) and (12) define a 
fermionic encoding that protects perfectly from the error, allowing the per- 
fect recovery of the information. The extension of the above procedure to 
the case of errors acting at different times {tj} is straightforward, by defining 
the sets {Fi(ti)} and imposing for each t% conditions (12). Suppose n is the 
number of distinct sites affected by the error £. There are 4n + 6 conditions 
for the AD parameters defining the q a 's. The minimal D for which a solution 
may be found is D = n + 2 , with 2 < D < N. Under this assumption, 
conditions (12), involving operators representing the error-action just on the 
decoding region, suffice to define fermionic modes globally unaffected by the 
error. 

In a practical sense, the procedure is worthwhile only if n <C N, since 
this confines the encoding and decoding procedures to the ends of the chain, 
preserving the central tenet of the state transfer protocol. We therefore 
take the condition n <C N as defining the concept of a low error rate in 
this scenario. It is important to emphasise, however, that this counting is 
in terms of the fermionic description of the error operators. Such counting 
is favourable for errors such as Z% or XiXi+i, but there are other local Pauli 
errors, such as Xi, which are necessarily described in terms of 0(i) fermionic 
modes, and therefore may not be included in the low rate condition. 

As previously mentioned, the typical assumption that the chain is ini- 
tialised in some global initial state \ipo) is not necessary [4, 19, 20], and so 
it is in our procedure. Indeed, there are two contexts in which our results 
apply equally. Either we can choose what the initial state should be, as we 
have so far assumed, or we are given a fixed initial state, in which case (10) 
gives two additional constraints on the encoding operators, which can be 
accommodated by our choice of D. The advantage of choosing a particular 
initial state is that it may reduce the effective number of errors that we have 
to correct for. The following section provides an illustration of this idea. 

4 Determining the Error 

Given a spin-chain affected by an (unknown) error £ acting at time t, con- 
ditions (12) indicate what knowledge of the error is needed to apply our 
encoding procedure. What experiment can one perform in order to gain 
this knowledge? One may, of course, use process tomography [ ] to deter- 
mine £ (tf — t), but this is extremely inefficient and would provide a lot of 
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redundant information. Our description of the error suggests indeed that 
there must be more efficient procedures, based on preparing the chain in a 
suitable set of states and then applying state tomography just on the decod- 
ing region to determine the error modes {F^}, with little regard for those 
acting on the complement region, {F^}- We shall now define a probing pro- 
cedure in the case = |0), motivating the existence of such procedures 
by using the most commonly assumed fixed initial state. 

In the case where the entire chain is initialised in the state = |0), 
it is only necessary to reconstruct a subset of the error modes. Indeed, 
observe that the encoding operators contain only creation operators, i.e., 

rjf = 0, Vi = 1, . . . D, hence in (12) the equations {q a ,F k i} = 0, V£,j are 

i 

automatically satisfied. Also, note that in £ the fermionic operators are 
ordered so that the annihilation operators act first. If |V>o) = |0)j there 
are no excitations on the complement region and therefore all error strings 
with any annihilation operators on the complement region just give a 
zero-contribution to the final state. In other words, from Eq. (11), for any 
i, \4> l s ) 7^ if and only if x G C, with C = {x : Xj = , V j > rrii + 1}. These 

x correspond to including no Ffc. Furthermore, Vx G C, since ql \0}® D 
belongs to the 1-excitation subspace, only the error modes containing no 
more that one annihilation operator F^ acting on the decoding region do 
not annihilate the encoding state. Namely, PxQa \ 7^ if and only if 
rii = mi or m = mi + 1. Let us define the sets S a = {i : rii — rrii = a}. We 
have just shown that to protect against the error it is sufficient to protect 

(i) 

against all the Fx such that i G S a , a = 0, 1, and x G C. So, it is sufficient 
that the encoding operators satisfy: 

{qLFv.} = 0, Vj = 1 . . . rij, Vi G So 

j 

{qai F k*} = 0,Vj = 1 . . .ni,Vi G 5i. 

j 

Moreover, Vz G Si, the latter reduces to just 

{qt,F kl } = 0,VieSi 

because this implies that pj q\ |0) 8Z) = , Vx G C, so imposing the condi- 
tions for j = 1 . . . rrii becomes unnecessary. 

In order to provide the information relevant to these conditions, we will 
probe the system with two different initial states, |0) = |0)®^ and |1) = 
\1)® D \0)® D . Note that both of these are prepared using the same fixed 
initial state outside the encoding region (in which we have the ability to 
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prepare any state), and this fixed state is the same one that will be used 
for the state transfer protocol. Consider first using |0). The only error 

(i) 

operators that do not annihilate the probing state are the P x such that 
i £ Sq and x £ C, which indeed do not include any annihilation operators. 
In order to reconstruct the information about the corresponding , we 
can therefore post-select on the decoding region being the the one-excitation 
subspace, and perform tomography to determine the corresponding density 
matrix. The span of states is described by the set Vq = {F k _ |0)® D }, Vj = 
1 . . . nii,Vi £ Sq. Hence, state tomography [2 I] applied just to the decoding 
region allows one to reconstruct SpanjVo}, which is sufficient information 
to impose {ql, F k i } = 0, Vj = 1 . . . rij, V? £ So- 

3 

When we use the initial state |1), the only non-zero contribution to the 

(i) 

final state is given by all the P x \/x £ C, i = 1 . . . s. Hence, the projection 

of this state on the D — l excitation subspace of the decoding region at time 

U) (i) 
tf includes the action of the Px with x £ C and i £ S Wx+ i, i.e., the P x 

where the one more operator acts than Fj operators. The span of these 

states hence includes those we are interested in, V\ = {-Pfc m . +1 \l)® D }ieSi j 
but also includes some others, V . This information is sufficient to impose 
the remaining conditions {qt,F k i } = 0, Vi £ Si, as desired, together with 
additional (unnecessary) conditions, which would protect from the error- 
components contributing to V . This redundancy is well tolerated, however, 
since we assume that the number n of the single error-modes is small com- 
pared to N. 

Overall, this procedure gives a method to determine the relevant errors 
using only operations on the encoding/decoding regions, and requiring a 
number of measurements that scales as 0(2 D ), and D ~ n <C N, thereby 
achieving a significant efficiency saving compared to standard process to- 
mography on the whole chain. 

5 Conclusions 

We have presented a protocol, incorporating encoding and decoding pro- 
cedures, to achieve perfect state transfer in the presence of low rates of 
systematic errors, whose repeatability allows one to learn about their struc- 
ture. This procedure may be thought of as an error correction optimally 
tuned to the error. Indeed, it ensures perfect recovery and, in addition, it 
has the appealing feature that if the number of sites affected by the error 
is small compared to the dimension of the whole chain, the encoding and 
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decoding operators involve just a small number of qubits, preserving the 
central feature of the ideal state transfer protocol. 

The low rate limit may be relaxed slightly, by extending the procedure 
to include errors with a support greater than the upper bound n, by se- 
lecting, via the proposed probing procedure, the first n operators Fk having 
the highest probability of acting on the decoding region and then encod- 
ing against them. This would minimise the probability of an unrecoverable 
error. 

Our formalism can also be applied in more general scenarios, where the 
action of the error is not instantaneous, since it is sufficient to describe the 
effect of any error just at the output time, tf. For instance, consider the 
case of a perturbed Hamiltonian, such as H = Hq + 5V, where Hq is the 
perfect state transfer Hamiltonian while 6V(t) is a perturbation which acts 
non-trivially only for a short time interval 5 at time t. In the interaction 
picture, the dynamics is determined by H (t) = U(t)5VU^ (t). Hence, 

r 5 

£ = U(t f - t) expH / H'(T)dT)U^(t f - t) |Vw) • 
J o 

For small 5, exp(— i H' (t)<1t) can be written just like the error (1), af- 
fecting, according to the Lieb- Robinson bound [ ] , only a small number of 
sites localised around the region where the perturbation acts. Therefore, the 
proposed encoding may equally be applied to recover perfect state transfer 
in the presence of this class of perturbations. 

In the future, it will be interesting to see how the error correcting ca- 
pabilities can be developed further, either by encoding on a single chain 
or across multiple chains, or by using different network topologies for com- 
munication. It will also be important to understand if errors such as local 
bit-flips can be corrected for. 
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A Construction of the decoding unitary 

To explicitly construct the decoding unitary, one orthogonalises the vectors 
belonging to each Z a via the Gram- Schmidt procedure: 

(iA _ y (^ 



The number z of such vectors satisfies z < Y^j=i 2™ J ( as some vectors in Z a 
may be linearly dependent). Then, Ud can be defined as: 

Ud = ^2 (\ w > ^ ® (* w l\ ' 

w,l a 

where the sum runs over all the orthogonalised vectors and {\w,h)} is any 
set of orthonormal vectors, defined over all the decoding region but qubit 
N. 

To see that the above is the desired unitary, one rewrites |V>f) in terms 
of the orthogonalised vectors and uses 

[6 Wtl \0i 3) ) = (vH 1 ^) = Vj,x,i,y, 
which holds because (6), (10) and (13) imply 

(l«|lW) = (0«|0W) , V^,Vy, Vi,j. 
Consequently, Ud applied to \^f), gives: 

\^ d ) = (t D ® U D ) |Vf> = m ® (a\0) N + ^1)^) , 
where |$) = X^e{o,i}™ 4 l^s) ® \ w ^)- Since 5 is trace-preserving, 

Th..JV-l(|^d> (^I) = |X> (Xl> 

i.e., the information has been perfectly transferred to the last qubit, as 
promised. 
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